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Abstract 

The ZEUS detector has been used to study dissociation of virtual photons in 
events with a leading proton, '~f*p — > Xp, in e~^p collisions at HERA. The data 
cover photon virtualities in two ranges, 0.03 < < 0.60 GeV^ and 2 < < 
100 GeV^, with Mx > 1-5 GeV, where Mx is the mass of the hadronic final state, 
X. Events were required to have a leading proton, detected in the ZEUS leading 
proton spectrometer, carrying at least 90% of the incoming proton energy. The 
cross section is presented as a function of t, the squared four-momentum transfer 
at the proton vertex, $, the azimuthal angle between the positron scattering 
plane and the proton scattering plane, and Q^. The data are presented in terms 
of the diffractive structure function, F^^^\ A next-to-leading-order QCD fit to 
the higher-Q^ data set and to previously published diffractive charm production 
data is presented. 
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1 Introduction 



In diffractive processes in hadron-hadron or photon-hadron collisions, the initial state 
particles undergo a "peripheral" collision, in which either the particles stay intact (elastic 
scattering), or they dissociate into low-mass states (diffractive dissociation). The scattered 
hadrons (or the low-mass states in the dissociative case) have energy equal, to within a few 
per cent, to that of the incoming hadron, and very small transverse momentum. Diffractive 
interactions can be parameterised in the framework of Regge phenomenology, where they 
are ascribed to the exchange of a trajectory with the vacuum quantum numbers, the 
Pomeron trajectory [1]. In the same framework, events in which the proton loses a more 
substantial fraction of its energy are ascribed to the exchange of subleading trajectories. 

Significant progress has recently been made in understanding diffraction in terms of QCD, 
notably by studying the diffractive dissociation of virtual photons in electron-proton or 
positron-proton collisions at HERA. In fact, diffraction has proven to be a tool to study 
QCD and the low- a; structure of the proton [2,3]. In the proton's rest frame, diffractive 
ep scattering, ep — > eXp, proceeds from the fluctuation of the virtual photon emitted 
by the electron (or by the positron) into a colour dipole, such as a quark- ant iquark pair 
or a quark-antiquark-gluon system. The dipole interacts hadronically with the proton 
via the exchange of an object with vacuum quantum numbers - a gluon pair, in leading- 
order QCD - and then dissociates into the hadronic state, X. The dipole has transverse 
dimensions which decrease as the photon virtuality, Q^, increases. It is thus possible to 
study diffractive interactions in a regime where one of the two interacting hadrons is so 
small that the strong interaction can be treated perturbatively. Alternatively, in a frame 
in which the proton is fast, the reaction can be seen as the deep inelastic scattering (DIS) 
of a pointlike virtual photon off the exchanged object. This gives access to the diffractive 
parton distribution functions (PDF), for which a QCD factorisation theorem has been 
proven [4-6]. Diffractive PDFs are defined as the proton PDFs probed when the vacuum 
quantum numbers are exchanged and the proton emerges intact from the interaction, 
suffering only a small energy loss. In the context of QCD, Pomeron exchange should then 
be understood as a synonym for exchange of partons from the proton with the vacuum 
quantum numbers. 

This paper presents new measurements of the reaction e^p — > e^Xp in the regions 0.03 < 
< 0.60 GeV^ (low-Q^ sample) and 2 < < iqo GeV^ (high-Q^ sample) ^ The 
measurements were made using the ZEUS detector at the HERA ep collider. The events 
were selected by requiring the detection in the ZEUS leading proton spectrometer (LPS) 
of a scattered proton, carrying a fraction Xi of the incoming proton momentum of at 

^ In the following, for simplicity, the symbol e will be used to denote both electrons and positrons. 
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least 0.9; such fast protons are referred to as leading. This xl range includes the so-called 
diffractive peak, the narrow peak in the cross section at xl ~ 1 ascribed to Pomeron 
exchange, as well as the transition to the lower xl region in which subleading, mesonic 
exchanges, notably Reggeons, are thought to dominate [1,3]. The measurement covers the 
region 0.075 < \t\ < 0.35 GeV^, where t is the square of the four-momentum transferred 
at the proton vertex. 

Sections 2-9 present the experimental set-up and the details of the analysis. Section 10 
gives the results. The t and $ dependences of the cross section are discussed first, where 
$ is the azimuthal angle between the positron and the proton scattering planes in the 7*p 
rest frame. The distribution of $ is sensitive to the interference between the amplitudes 
for scattering of longitudinally and transversely polarised photons, and thus to the helicity 
structure of the interaction. The dependence of the photon-proton differential cross- 
section da'^'P^'^P/dMx is then studied for different values of the photon-proton centre-of- 
mass energy, W. The data are also discussed in terms of the diffractive structure function, 
^i)(3) rpj^^ dependence of F^^^^ on xjp is studied, where xjp is the fraction of the proton 
momentum carried by the object (the Pomeron or the Reggeon, in the Regge framework) 
exchanged between the virtual photon and the proton, xp ~ 1 — x^. The and f3 
dependences of F^^^^ for different values of xjp are investigated, and the behaviour of 
^2^^^^ is compared to that of the inclusive proton structure function, F2. The variable (3 is 
the Bjorken variable defined with respect to the four-momentum of the exchanged object. 
The results are compared to theoretical predictions based on the colour-dipole approach 
outlined above. Finally, a next-to-leading-order (NLO) QCD fit to the higher-Q^ data is 
presented. 

The present data correspond to an integrated luminosity about a factor four larger than 
that of the previous ZEUS-LPS analysis [7]. The low-Q^ results have a wider and W 
coverage than that studied so far with the LPS or other methods [7,8]. The $ distribution 
is investigated for the first time. The cross section is measured up to Mx values of 
40 GeV, so far unexplored, and ^2^^^^ is presented up to xp values of 0.07, thus covering 
the diffractive-peak region (xjp < 0.02) and the transition to the non-diffractive region 
that is dominated by subleading exchanges. 

It is also possible to select diffractive events without detecting the scattered proton. In a 
previous paper [9], a method based on features of the shape of the mass spectrum of the 
hadronic final-state X {Mx method) was applied. A discussion of the two approaches is 
presented, along with a comparison of the corresponding results. 
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2 Kinematics and cross sections 



Figure 1 shows a schematic diagram of the process ep —> eXp. The kinematics of this 
reaction is described by the variables: 

• = —q^ = —{k — k')^, the negative four-momentum squared of the virtual photon, 
where k {k') is the four-momentum of the incident (scattered) positron; 

• W"^ = (q + P)"^, the squared centre-of-mass energy of the photon-proton system, where 
P is the four-momentum of the incident proton; 

• X = I {2P ■ q), the fraction of the proton momentum carried by the quark struck by 
the virtual photon in the infinite momentum frame (the Bjorken variable); 

• M\ = {q+P-P'f, the squared mass of the system X, where P' is the four-momentum 
of the scattered proton; 

• t = {P — P'Yi the squared four- momentum transfer at the proton vertex; 

• $, the angle between the positron scattering plane and the proton scattering plane in 
the 7*p centre-of-mass frame. 

The variables Q^, W and x are related by x = Q'^ /{Q'^ + — M^), where Mp is the 
proton mass. 

The differential cross sections for the reactions ep eXp and 7*p Xp are related by 



1 - 



l + {l-y)n + RD 



dlnW^dQ^dMxdt 27r dMxdt 

a l + (l-y)^rfaT>-^P 

27r Q2 dMxdt ' ^ ' 

where a is the fine structure constant and y = {P ■ q)/{P ■ k) is the fraction of the positron 
energy transferred to the proton in its rest frame. The quantity = aj^ P~*^P/a'^ p^Xp 
the ratio of the cross sections for longitudinally and transversely polarised virtual photons. 
In the region covered by the present data, the term within the square brackets is taken 
to be unity since R^ is expected to be small [10]. 

The kinematics of the reaction ep —>■ eXp can also be described by Q^, t and $, in 
conjunction with the two dimensionless variables xjp and (3 introduced in Section 1 given 
by 

_ {P-P')-q _ Q^ + Ml-t 
P-q ~Q^ + W^-Ml' 

3 (3) 
^ 2{P-P')-q Q^ + Ml-t ^' 

The quantities xjp and (3 are related to x by = x. 
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The cross section for the reaction ep — > eXp can be expressed in terms of the structure 
function F^^^\ which is defined by the equation 



d(5dQ'^dx]pdt 



y 



2{1 + R^) 



V + — 
^ 2 



F,''^'\(3,Q'\xjp,t). 



(4) 



The structure function F^^^\j3,Q'^,xip) is obtained by integrating F^^'^' over t: 



D(4) 



F,''^'\p,Q',xjp) 



,D(4) 



{p,Q\x]p,t)dt. 



For the results presented in this paper, the integration was performed in the range < 
\t\ < 1 GeV^. 



For unpolarised positrons and protons, the cross section can also be decomposed as 



oc 4*^^"^^ 



ecr 



e a 



LT 



COS $ — ea, 



"f'p—tXp 
TT 



COS 2$, 



(5) 



where a^^^^^ is the interference term between the amplitudes for longitudinal and trans- 
verse photons and (Jrprp IS the interference term between the amplitudes for the two 
transverse polarisations. The polarisation parameter e is defined as e = 2(1 — y)/[l + (1 — 

y?]- 



3 Selection of diffraction at HERA 



The kinematics of diffractive scattering, 7*p Xp, implies that three features should be 
present in the final state: 

1. the proton suffers only a small perturbation and emerges from the interaction carrying 
a large fraction, xl, of the incoming proton momentum. Diffractive events appear as 
a peak at ~ 1, the diffractive peak, which at HERA approximately covers the 
region 0.98 < Xl < I- The absolute value of the four-momentum transfer squared t is 
typically much smaller than 1 GeV^, with {\t\) ~ 0.15 GeV^; 

2. conservation of momentum implies that any other produced system (X) must have 
a small mass (Mx) with respect to the photon-proton centre-of-mass energy (since 
1-xl>M],/W''); 

3. the difference in rapidity between the outgoing proton and the system X is Ar] ^ 
In (1/xjp) [3]. This, combined with the peaking of the cross section at small values of 
xjp, leads to a large separation in rapidity between the outgoing proton and any other 
hadronic activity in the event. 
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There are two basic ways to select inclusive diffractive events. The first is the proton- 
tagging method (exploiting the first signature above), used in the present study and in 
earlier ones [7,8,11-13]. The second exploits the different characteristics of the system X 
in diffractive and non- diffractive events: 

• in non- diffractive DIS, both the hadronic system associated with the struck quark, 
which is largely measured in the detector, and that of the proton remnant, which 
largely escapes down the beam-pipe, are coloured states. In this case, the distribu- 
tion of the final-state particles is governed by conventional QCD fragmentation and 
particles are expected to be emitted uniformly in rapidity along the '~f*-p axis. This 
leads to a suppression of rapidity gaps as well as a suppression of small masses of the 
hadronic system observed in the detector; 

• in contrast, small masses of the system X and large rapidity gaps are signatures of 
diffractive processes (the second and third signatures above). At HERA, diffractive 
analyses based on the hadronic methods have been made with event selections based 
both on the presence of large rapidity gaps (rapidity-gap method, see e.g. [2, 14] and 
references therein) and small masses of the system X {Mx method) [7,9]. 

The two basic approaches for the selection of diffractive events, the proton-tagging method 
and the hadronic methods, are complementary: 

• in the hadronic methods, high Mx values are not accessible since the non-diffractive 
background grows with Mx- Also, the estimation of the non-diffractive background 
relies on models of fragmentation. Furthermore, the measured cross section includes 
a contribution from proton-dissociative events, ep — eXN. In these events, the mass 
of the forward system (A^) enters as another variable, and the observed particles must 
be assigned either to the system A^ or X. The number of events in which no particle 
from A^ is observed must be estimated from a Monte Carlo simulation. While these 
limitations add to the systematic uncertainty of the hadronic methods, the statistical 
precision of the results tends to be good due to the high acceptance of the central 
detector. Also, the acceptance is not limited in t, although no measurement of t is 
possible; 

• conversely, samples selected in the proton-tagging method have little or no background 
from proton-dissociative events or from non-diffractive DIS. They also allow a direct 
measurement of the variables t, $ and xp (at large values of xp), and give access 
to higher values of Mx- The statistical precision, however, is poorer than for the 
results obtained using hadronic methods due to the small acceptance of the LPS - 
approximately 2% in the diffractive peak region. 

Section 10.5 presents a comparison between the results obtained with the LPS and Mx 
methods. The results are also compared to measurements made by the HI collabora- 
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tion [14] in which diffraction was selected with the rapidity-gap method. 



4 Experimental set-up 

The measurements were carried out at the HERA colhder in 1997 using the ZEUS detector. 
At that time, HERA operated with 820 GeV protons and 27.5 GeV positrons. The data 
used in this analysis correspond to integrated luminosities of 3.60 ± 0.06 pb~^ (low-Q^ 
sample) and 12.8 ± 0.2 pb~^ (high-Q^ sample). 

A detailed description of the ZEUS detector can be found elsewhere [15]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [16], which operates 
in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD consists 
of 72 cylindrical drift chamber layers, organised in nine superlayers covering the polar- 
angle^ region 15° < 9 < 164°. The transverse-momentum resolution for full-length tracks 
is (t{pt)/pt = 0.0058pT © 0.0065 © 0.0014/pr, with pr in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [17] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. The CAL 
energy resolutions, as measured under test beam conditions, are a(E)/E = 0.18/ ^/E for 
electrons and a(E)/E = 0.35/^/E for hadrons {E in GeV). 

Low-Q^ events were selected by identifying and measuring the scattered positron in the 
beam-pipe calorimeter (BPC) [18] and beam-pipe tracker (BPT) [19]. The BPC was a 
tungsten-scintillator sampling calorimeter, located 3 m from the interaction point and 
covered positron scattering angles relative to the incident direction of 15 to 34 mrad. The 
BPT was a silicon-microstrip tracking device situated immediately in front the BPC. In 
1997, it was equipped with two detector planes to measure the X coordinate. 

For the high-Q^ sample, the impact point of the scattered positron was determined with 
the small-angle rear tracking detector (SRTD) [20] or the CAL. The SRTD is attached 
to the front face of the RCAL and consists of two planes of scintillator strips, 1 cm wide 
and 0.5 cm thick, arranged in orthogonal orientations and read out via optical fibres 
and photomultiplier tubes. It covers a region 68 x 68 cm^ in X and Y, excluding a 
10 X 20 cm^ hole at the centre for the beam-pipe. The corresponding angular coverage is 
between 4° and 18° around the beam-pipe. Ambiguities in SRTD hits were resolved with 
the help of the hadron-electron separator (HES) [21], which consists of a layer of 10 000, 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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2.89 X 3.05 cm^ silicon-pad detectors inserted in the CAL at a depth of 3.3 radiation 
lengths. 

The LPS [22] detected positively charged particles scattered at small angles and carrying a 
substantial fraction, x^, of the incoming proton momentum; these particles remain in the 
beam-pipe and their trajectory was measured by a system of silicon microstrip detectors 
that could be inserted very close (typically a few mm) to the proton beam. The detectors 
were grouped in six stations, SI to S6, placed along the beam line in the direction of 
the proton beam, between 23.8 m and 90.0 m from the interaction point. The particle 
deflections induced by the magnets of the proton beam-line allowed a momentum analysis 
of the scattered proton. For the present measurements, only stations S4, S5 and S6 were 
used. The resolutions were about 0.5% on the longitudinal momentum and about 5 MeV 
on the transverse momentum. The effective transverse-momentum resolution is dominated 
by the intrinsic transverse- momentum spread of the proton beam at the interaction point, 
which is about 40 MeV in the horizontal plane and about 90 MeV in the vertical plane. 
The LPS acceptance is approximately 2% and x^-independent in the diffractive-peak 
region, 0.98 < xl < 1; it increases smoothly to about 10% as xl decreases to 0.9. 

The luminosity was measured from the rate of the bremsstrahlung process ep e'-fp. The 
photon was measured in a lead-scintillator calorimeter [23] placed in the HERA tunnel at 
Z = -107 m. 



5 Reconstruction of the kinematic variables 

In the low-Q^ analysis, 0.03 < < 0.60 GeV^, the scattered positron was measured in the 
BPC/BPT. The energy and angle of the scattered positron were used ("electron method") 
to determine Q'^ and W. For the high-Q^ data (2 < < 100 GeV^) the identification 
of the scattered positron was based on a neural network [24] which uses information from 
the CAL. The variables W and were reconstructed using a combination of the electron 
method and the double angle method [25]. 

The longitudinal (pz) and transverse (pxiPy) momenta of the scattered proton were 
measured with the LPS. The fractional energy of the outgoing proton, xl, was defined as 
Xl = Pz/Ep, where Ep is the incoming proton energy. The variable t is given by 

t = -?i-ii^iV$, (6) 

XL Xl ^ 

where px is the transverse momentum of the proton with respect to the incoming beam 
direction. The t resolution is approximately (T{\t\) = 0.14 GeVA/]t| {t in GeV^) and is 
dominated by the angular spread of the beam. The proton and the positron momenta 
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were used to determine $, the azimuthal angle between the positron and proton scattering 
planes in the 7*p frame. The resolution on $ is approximately 0.2 rad. 

The four-momentum of the system X was determined from calorimeter and tracking 
information. The energy deposits in the CAL and the track momenta measured in the 
CTD were combined in energy flow objects (EFOs) [9, 26] to obtain the best momentum 
resolution. The mass Mx was evaluated as 

mIj,fo=(J2^) -{J2p^^) -{T.py^) -(JIp^^^ 

where the sums run over all EFOs not assigned to the scattered positron. The mass Mx 
can also be determined from the outgoing proton momentum as reconstructed in the LPS: 

The best resolution on Mx is obtained with Mx,efo when Mx is small and with Mx,lps 
for large values of Mx, so Mx was reconstructed as 

M|. = wefoM^ 

EFO + "^LPsMx rpg, (7) 

where the weights wefo and wlps are inversely proportional to the appropriate resolu- 
tions, and wefo + wlps = 1- The resulting resolution is a{Mx)/Mx = 0.35/a/Mx + 0.08, 
with Mx in GeV. 

The variables xjp and /5 were obtained from Eqs. (2) and (3), using the value of Mx 
from Eq. (7). For the determination of xp, this procedure is equivalent, at large Mx, to 
evaluating xjp as xjp = 1 — xl- 

The variable y was reconstructed as = X] *"^'~2g ' where the sum is over all EFOs 
not associated to the scattered positron and is the energy of the incident positron 
("Jacquet-Blondel method" [27]). 

6 Event selection 

The data used for the analysis were selected at the trigger level by requiring the presence 
of a scattered positron in the BPC or CAL and a scattered proton in the LPS. In the 
offline selection the following cuts were imposed, closely following those used in the F2 
analyses at low Q"^ [19] and high Q"^ [28]: 

• the energy of the scattered positron, if measured in the BPC, was required to be 
between 3 and 7 GeV for the events with W > 260 GeV and greater than 7 GeV for 
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W < 260 GeV, reflecting the trigger selection. If measured in the CAL, the scattered 
positron energy was required to be larger than 10 GeV. 

The position of the scattered positron was required to be within the fiducial regions 
of the BPC or the CAL. In addition, in case the positron was found in the BPC, the 
impact position at the BPC front face, as extrapolated from the BPT measurement, 
was required to match with that of the BPC shower. Furthermore, in order to identify 
electromagnetic showers and to reject hadrons, the transverse size (energy weighted 
root mean square) of the shower in the BPC was required to be less than 0.8 cm; 

• the requirements 30 < {E — pz) < 65 GeV and 40 < {E — pz) < 65 GeV were 
imposed for the low- and high-Q^ samples, respectively, where E — pz = Yl i^i ~ Vzi) 
and the summation runs over the energies and longitudinal momenta of the final-state 
positron and all EFOs. This cut reduces the size of the QED radiative corrections and 
the photoproduction background, i.e. the ~ events where the scattered positron 
escapes undetected in the rear beam-hole; 

• in order to limit event migrations from low the variable yjB was required to be 
greater than 0.06; 

• the Z coordinate of the interaction vertex was required to be in the range —90 < 
Z < 90 cm for the low-Q^ sample and —50 < Z < 50 cm for the high-Q^ sample. 
Events without a vertex reconstructed using BPT tracks were discarded in the low-Q^ 
analysis. In the high-Q^ sample, events without a measured vertex were assigned to 
the nominal interaction point. 

The following requirements were used to select the scattered proton measured in the LPS: 

• the candidate proton was tracked along the beam line and was rejected if, at any 
point, the distance of approach to the beam pipe was less than 0.2 cm for px < 0, or 
less than 0.3 cm for px > 0. This cut reduces the sensitivity of the acceptance to the 
uncertainty in the position of the beam-pipe apertures; 

• the variable t was required to be in the range 0.075 < \t\ < 0.35 GeV^. This cut 
eliminated regions where the LPS acceptance was small or rapidly changing; 

• beam-halo background is caused by scattered protons, with energy close to that of 
the beam, originating from the interaction of a beam proton with the residual gas 
in the beam-pipe or with the beam collimators. A beam-halo proton may overlap 
with a standard non- diffract ive DIS event. In this case, the proton measured in the 
LPS is uncorrelated with the activity in the central ZEUS detector. This background 
was suppressed by the requirement that the sum of the energy and the longitudinal 
component of the total momentum measured in the CAL, the BPC and the LPS 
be less than the kinematic limit of twice the incoming proton energy: E + pz = 
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{E + pz)cAL + {E + Pz)bpc + 2p^^^ < 1655 GeV. This cut takes into account the 
resolution of the measurement of p^^. The residual beam-halo background and its 
subtraction are discussed in Section 7. 

The low-Q^ analysis was further limited to the kinematic region 0.03 < < 0.6 GeV^, 
63 < W < 280 GeV, Mx > 1.5 GeV and xl > 0.9. The average value for this 
sample is 0.23 GeV^. The high-Q^ analysis was restricted to 2 < < 100 GeV^ 
25 <W < 240 GeV, Mx > 1.5 GeV and xp < 0.07; the average value is 10.5 GeV^. 
These selections yielded 334 events in the low-Q^ sample and 5945 events in the high-Q^ 
sample. 

7 Monte Carlo simulation and acceptance corrections 

Monte Carlo simulations were used to correct the data for acceptance and detector effects. 
In the low-Q^ analysis, events of the type ep — eXp were simulated with the generator 
EPSOFT 2.0 [29,30], based on the triple- Regge formalism [1,3], in which the cross section 
can be expressed in terms of three trajectories. If all the trajectories are Pomerons 
[IP IP IP), the cross-section da'^'^'^^^/dMjr is approximately proportional to l/Mj.. If 
one of the trajectories is a Reggeon (IPIPIR), the cross-section d(j"'*'P~*-^^ / dM'j^ falls as 
~ I /Ml. 

In the high-Q^ analysis, the reaction ep — > eXp was modelled with RAPGAP 2.08/06 [31], 
which is based on the model of Ingelman and Schlein [32]. In RAPGAP, the structure 
function -Fg^^^"* is expressed as the sum of Pomeron and Reggeon contributions: 

F^'^^ {xjp, t, /3, = fjp{xjp, t)Ffi/3, Q') + fnixjp, t)F^{P, Q'). (8) 

The Pomeron and Reggeon fluxes, fip,iR{x]p,t), were parameterised [1] as 

,iP,K, 

Xjp 

with linear trajectories ajp^nit) = ajp^mi^O) + Oi'jpjpit, and the values of the parameters were 
taken from hadron-hadron data. The Pomeron structure function F^{j3, Q"^) was taken 
from the HI measurements [14] (fit 2). The Reggeon trajectory includes the p, the u, the 
/2 and 02 mesons. Their structure functions are unknown and were approximated with 
that of the pion [33]. The assumption that F^^'^^ can be expressed as the product of a 
flux, depending only on xjp and t, and the structure function of a particle-like object (see 
Eq. (8)) is known as the "Regge factorisation" hypothesis. It gives a fair description of 
the data, although it has no justification in QCD, where only the concept of diffractive 
PDFs, which are functions of (xjp, t, /3, Q^), has a firm basis. 
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The Monte Carlo generator DIFFVM 1.0 [34] was used to simulate the double-dissociative 
reaction, ep — > eXN, where the proton diffractively dissociates into the state N. 

Initial- and final-state QED radiation were simulated by using EPSOFT and RAPGAP 
in conjunction with HERACLES 4.6 [35]. 

The generated EPSOFT and RAPGAP events were reweighted in xl and t so that the 
measured distributions were well described. All generated events were passed through the 
standard ZEUS detector simulation, based on the GEANT program [36], and the trigger 
simulation package. A comparison of data and MC simulations is presented in Figs. 2 
and 3 for the variables xl, t, Q^, W, Mx, and xjp. In Figs. 2c and 2d, no LPS requirement 
was imposed, so as to reduce the statistical fiuctuations. The simulations reproduce the 
data satisfactorily. The diffractive peak is evident in Figs. 2a and 3a. 

In the low-Q^ analysis, the measured number of events was corrected for acceptance bin- 
by-bin. Bin-centring corrections were applied in W, xjp and Mx assuming u'^'p^^p oc 
W^'^'^'''P-^\ F^''^^ oc l/xp'^~\ and da^'P^^P/dM]. oc 1/M^^^"^"^\ where ajp is the t- 
averaged value of the Pomeron trajectory as obtained from the present measurement 
(Section 10.6). In the high-Q^ analysis, the cross section for the dissociation of virtual 
photons at a given point within a bin was obtained from the ratio of the measured number 
of events to the number of events in that bin predicted from the MC simulation, multiplied 
by the 7*p — > Xp cross section calculated by the Monte Carlo generator. Both the 
acceptance and the bin-centring corrections were thus taken from the MC simulation. 

The cross section was directly measured only in a limited t region and extrapolated to 
< |t| < 1 GeV^ using the t dependence assumed in the Monte Carlo generator, which was 
reweighted to the measured value of b (see Section 10.1). The effect of the extrapolation is 
to increase the cross section by a factor of about two (to within 3%); this factor is largely 
independent of the measured kinematic variables. In the region covered by the present 
measurements, the extrapolation is performed assuming an exponential dependence on 
t, da'^P^'^^P / dt oc exp(6t), with b ~ 7.9 GeV~^ at low X]p. Data from elastic and p- 
dissociative pp and pp scattering indicate that the t distribution is better described by 
the function exp {bt + ct^). For example, fits to the pp data at ^/s = 546 GeV [37] yield 
c = 2.3 ± 0.1 GeV~'^. The present data are only weakly sensitive to the value of c. In the 
extrapolation to the range < |t| < 1 GeV^, c was taken to be zero, but was varied by 
up to 4 GeV""^, yielding changes in the extrapolated cross section of up to 7%. This was 
included in the normalisation uncertainty discussed in Section 9. 

The results presented in this paper were corrected to the Born level and to the following 
kinematic region: 

• low-g2 data: 0.03 < Q"^ < 0.6 GeV^, 63 < W < 280 GeV, Mx > 1.5 GeV, xl > 0.9 
and < |t| < 1 GeV^; 
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• high-g2 data: 2 < < iqo GeV^, 2b < W < 240 GeV, Mx > 1.5 GeV, xl > 0.9 
and < Itl < 1 GeV^. 



8 Backgrounds 

The main background contribution is given by proton beam-halo events. In such events, 
the proton detected in the LPS is not correlated with the measurements in the central 
detector. To estimate the size of this background, the variable E + pz (see Section 6) was 
used. For a signal event, this quantity should be equal to twice the initial proton energy, 
1640 GeV, whereas for a beam-halo event it can exceed this value. 

The E + pz spectrum for the beam-halo events was constructed as a random combination 
of generic DIS events (without the requirement of a track in the LPS) and a beam-halo 
track measured in the LPS. The resulting distribution, shown in Fig. 4 as the hatched 
histogram, was normalised to the data for E +pz > 1685 GeV, which contain beam-halo 
events only (see Fig. 4). The background remaining after the E + pz < 1655 GeV cut 
averaged (10.2 ± 0.7 (stat.))% and was a function of xp and t. All results presented in 
this paper were corrected for this background. 

The contribution from proton-dissociative events, ep — >■ eXN, was studied with the Monte 
Carlo generator DIFFVM for the DIS sample and was found to be less than 4% in all 
bins. This contribution was neglected. 

9 Systematic uncertainties 

The systematic uncertainties were calculated by varying the cuts and by modifying the 
analysis procedure. The variations of the cuts were typically commensurate with the 
resolutions of the relevant variables, and in general the changes were similar to those 
made in earlier analyses [11,12,19,22]. 

The following systematic checks were performed ^: 
• to evaluate the uncertainties due to the measurement of the scattered positron: 

— low-Q^ analysis: the checks performed in a previous publication [19] were repeated, 
and consistent results obtained. Since the present data are a subsample of those 

^ The corresponding average effect on the cross section in the measured bins is indicated using the 
notation given a systematic check which produces an increase of the cross section in some bins 
and a decrease in some other bins, a is the average increase and b is the average decrease. 
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used earlier, the systematic uncertainties found previously [19] were used (typically 
smaller than ±1.5%); 

— high-Q^ analysis: the fiducial region for the impact position of the positron was 
modified (Hg^o); the minimum energy of the positron was increased to 12 GeV 

to evaluate the uncertainties due to the reconstruction of the final-state proton, the 
checks described below were performed. Consistent results were obtained for the low- 
and the high-Q^ samples; because of the larger statistical fiuctuations of the low-Q^ 
sample, the uncertainties quoted are those determined from the high-Q^ events; 

— the cut on the minimum distance of approach to the beam-pipe was increased by 
0.03 cm (ttiVo); the t range was enlarged to 0.07 < \t\ < 0.4 GeV^ the 
amount of the subtracted beam-halo background was varied by ±20% (^i 4%); 

sensitivity to the other selection cuts: 

— low-Q^ analysis: the checks performed in a previous publication [19] were repeated 
and consistent results found, notably for the sensitivity to the selections on yjs, 
E — pz and the Z coordinate of the vertex. The systematic uncertainties deter- 
mined in [19] were used (typically smaller than ±1.5%); 

— high-Q^ analysis: the minimum value of E — pz was raised to 45 GeV (^21%)) 
the cut on the Z coordinate of the vertex was restricted to —40 < Z < 40 cm 
(li'3%); the systematic error due to the uncertainty in the absolute calorimeter 
energy calibration was estimated by changing the energy scale by ±2% (li'|%); 
the minimum value of Mx was decreased to 1.3 GeV (1l4%); 

sensitivity to the Monte Carlo simulations: 

— low-Q^ analysis: no reweighting of the xl distribution was applied in the Monte 
Carlo simulation (Igs^); the value of the t slope was changed by ±1.5 GeV"^ 

— high-Q^ analysis: the xjp distribution was reweighted by a factor (1/xjp)'^', with k 
varying between —0.05 and ±0.05 (]'^}'i%); the value of the t slope was changed 
by ±1.5 GeV~^ (-21%)^ the $ distribution was reweighted by a factor (1 ± 
k cos $), with k varying between —0.15 and ±0.15 (li'2%); the intrinsic transverse- 
momentum spread of the proton beam at the interaction point was increased by 
10 MeV in the horizontal plane and 20 MeV in the vertical plane (li'^'l'^o); 

the t slope was determined with an alternative method [22] based on expressing the 
measured t distribution as a convolution of an exponential distribution, da/dt oc e"^'*', 
and a two-dimensional Gaussian distribution representing the transverse momentum 
distribution of the beam. This led to changes in the value of the t slope by up to ±4%. 
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The total systematic uncertainty for each bin was determined by quadratically adding the 
individual contributions. The quoted uncertainties do not include an overall normalisation 
uncertainty of ±10% which originates mostly from the uncertainty of the simulation of the 
proton-beam optics; this uncertainty is largely independent of the kinematic variables and 
was therefore taken as a normalisation uncertainty. The ±10% normalisation uncertainty 
also includes the uncertainty on the integrated luminosity (±1.6%). In addition, a ±7% 
uncertainty is present in the cross section and structure function results, except those for 
da'^P^'^^P /dt, due to the extrapolation from the measured t region to < |t| < 1 GeV^ 
(see Section 7). Thus the overall normalisation uncertainty is l}o^. 



10 Results 

The results in this section are presented as follows. The differential cross-section da'^^^'^^^f dt 
in the region 0.075 < \t\ < 0.35 GeV^ is presented first. The data are then integrated over 
t and extrapolated to the range < |t| < 1 GeV^, as discussed in Section 7. The resulting 
cross sections are presented as a function of $ in Section 10.2, where the sensitivity of 
the present data to the helicity structure of the reaction ep — * eXp is discussed. The 
dependence of the cross section on Q'^ is presented in Section 10.3, and is compared with 
that of the total photon-proton cross section. 

In Section 10.4 the data are presented in terms of the diffractive structure function, F^^^\ 
This allows an interpretation based on the diffractive PDFs of the proton. The present 
data in the diffractive-peak region are compared to those obtained with the Mx method 
by ZEUS and with the rapidity-gap method by HI in Section 10.5. 

In Section 10.6, the xjp dependence of ^2^''^'* is used to extract the intercept of the Pomeron 
trajectory, ajp(O), the quantity that, in Regge phenomenology, determines the energy 
dependence of the total hadron-hadron cross section [1]. It is interesting to see if the value 
of «jp(0) in ep diffractive scattering at high is larger than that measured in hadron- 
hadron collisions, as expected in the framework introduced earlier, in which diffraction is 
due to the exchange of partons from the proton. In Section 10.7 the results are compared to 
some perturbative QCD (pQCD) models based on the dipole picture outlined in Section 1. 
Finally, an NLO QCD fit was performed. 

The results of this paper extend up to xjp ~ 0.07. In the following, the data for the 
diffractive-peak region are often contrasted with those at high xjp. For this purpose, the 
value X]p = 0.01 is chosen as the transition between the high- and low-Xip bins, such 
that the low-xjp bins are dominated by diffractive-peak events. This choice is somewhat 
restrictive, since the diffractive peak extends to Xip ^ 0.02, see Figs. 2a and 3a. In the 
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region xjp < 0.01, the contribution from non-Pomeron exchanges is less than 10%. The 
average value of xjp is 0.003 for xip < 0.01 and 0.043 for xjp > 0.01. 

Preliminary versions of the present results, along with details of the analysis, have been 
presented earlier [38-40]. 

10.1 t dependence 

Figure 5a presents the differential cross-section da'^^~*'^^^ / dt in the kinematic range 2 < 
< 100 GeV\ Mx > 1.5 GeV, xjp < 0.01. For this sample, (Q^) = 8.4 GeV^ and 
{(3) = 0.32. The value of the slope parameter, b, obtained from the fit with the function 
da^P^^^P/dt oc e-*l*l in the range 0.075 < |t| < 0.35 GeV^ is 

6 = 7.9 ± 0.5(stat.)+°|(syst.) GeV~l 

This agrees with and improves on the previous measurement of the diffractive slope of 
b = 7.2± l.l(stat.)^[!](syst.) GeV-^ [n]. 

The t distribution was studied in two bins, 2 < < 7 GeV^ and 7 < < 100 GeV^ 
The fitted values of 6 are 7.7±0.7(stat.)+°:?(syst.) GeV-^ and 8.0±0.8(stat.)l[J|(syst.) GeV~2, 
respectively. These results, presented in Fig. 5b together with the previous ZEUS mea- 
surements in photoproduction [12] and in DIS [11], show that b is independent of 
within the errors. This behaviour is expected in a QCD-based model [41]. 

However, the value of fe decreases with X]p, as shown in Fig. 5c (see also Table 1). Figure 5d 
shows the values of 6^2^ in bins of xjp, calculated from the fit to da'^P~*'^^P/dp'^, so that 
they can be compared with those of a previous ZEUS publication [8]. The relationship 
between t and p'^ is given in Eq. (6). In the region of the present data, the values of b and 
6p2 differ by less than their uncertainties. The present results are consistent with those of 
the previous ZEUS publication [8]. The slope reaches a minimum value for xjp ~ 0.05 
and then rises to 6p2 ^ 7 GeV^^ for higher xp values. 

In the dipole model, in which the virtual photon fluctuates into qq or qqg systems, the 
qqg contribution dominates for /5 < 0.2-0.3 [3]. Different t dependences for the qq and qqg 
regions are expected [41], with slopes higher by up to 3 GeV~^ at high (3 in the diffractive 
peak. In addition, in the qqg region, the slope b is expected to decrease logarithmically 
with X]p: b = b^ — 2a'jp Inxjp, the so-called shrinkage of the diffractive cone [3]. A decrease 
is observed in the data, and so is a significant (3 dependence. However, these dependences 
are visible only over an xjp range, xjp < 0.07, that goes beyond the diffractive-peak region 
(Fig. 5c), making a comparison with the predictions inconclusive. At higher xjp values, a 
rise of b with xjp is predicted by a Regge-based model [42], though at a rate smaller than 
that observed in the data. 
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10.2 Azimuthal asymmetry 



The azimuthal angle $ between the positron and proton scattering planes is sensitive 
to the helicity structure of the reaction ep —>■ eXp, as shown explicitly in Eq. (5). The 
analysis of the azimuthal distribution is limited to the high-Q^ data, since the statistics 
of the \ow-Q^ sample is too small. For this part of the analysis, a radial cut of 18 cm was 
imposed on the impact point of the scattered electron at the RCAL surface, along with 
the restriction > 4 GeV^. This reduces the $ dependence of the acceptance. 

Figures 6a and 6b show the $ distribution for the two ranges 0.00025 < xjp < 0.01 and 
0.01 < xjp < 0.07. Here again, xjp = 0.01 was chosen as the transition between the high- 
and low-Xjp bins such that the low-xjp bin is dominated by diffractive-peak events. The 
distributions were fitted to the form 

; oc 1 + Alt cos $ + Att cos 2$, 

where Alt and Att are proportional to a^[r^~''^^ and cxj/^^^, respectively. The values 
of the azimuthal asymmetries obtained in the fit are 

Alt = 0.009 ±0.073(stat.)l°:°^^(syst.), 

Att = 0.005 ± 0.074(stat.)l[j:|j^^(syst.) 

and 

Att = 0.007 ±0.048(stat.)l°;°^?(syst.), 

Att = 0.019 ± 0.046(stat.)l[J;J]|(syst.) 

for the ranges 0.00025 < Xp < 0.01 and 0.01 < Xjp < 0.07, respectively. The interference 
terms between the longitudinal and transverse amplitudes and between the two trans- 
verse amplitudes thus appear to be small in the measured kinematic range, both in the 
diffractive-peak region and at higher xjp values, suggesting that the helicity structure of 
the reaction ep — > eXp is similar for Pomeron and Reggeon exchanges. 

Figure 6c presents Alt as a function of xjp and Figs. 6d-f present Alt as a function of 
(3, t and for xjp < 0.01. The asymmetry is consistent with zero in all measured bins. 
The results are summarised in Table 2. 

The measured value of Alt can be compared with the results obtained in exclusive elec- 
troproduction of p° mesons, ep ep^p, in which the hadronic final state, X, consists of 
a p° meson only. In this case, Alt = -x/2e(l + e)- (r^o + S^ii) = -0.262 ± 0.038(stat.) ± 
0.068(syst.), where Tqq and rf^ are two of the spin-density matrix elements [43]. The 
present data show that the asymmetry is smaller for inclusive scattering than for exclusive 
p° electroproduction. 
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There are numerous pQCD-based predictions for the behaviour of A^t [44-47] in the 
diffractive peak, mostly concerning the high-/? region {/3 > 0.6-0.9), where the asymmetry 
is expected to be largest, reflecting the large expected value of this region was 

not accessible due to limited statistics. In all calculations, back-to-back configurations, 
i.e. Alt < 0, are favoured; the asymmetry is expected to be close to zero at low /5, in 
agreement with the present data. 

The measurement of the $ dependence can, in principle, be used to constrain the cross sec- 
tion of longitudinally polarised photons [45-47] , a quantity notoriously difficult to extract 
unless data at different centre-of-mass energies are available. No experimental results on 
cr2 ^^^P exist so far. The asymmetry Alt can be related to [47]; however, only for 
/5 > 0.8-0.9, beyond the region covered by the present data, is the determination of 
model-independent. More general limits can be obtained for a1 ^^"^^; using hermiticity 
and parity conservation, the following relations are found [45,46]: 



r) / 7 p — ^^P\ Z ) 



-/*p-^Xp-^2 „ _ „1*P^Xp 

'LT 

1*P 

ClQ — CItT 

where ctq = (Tt ^^^^ -\- eaj^ The larger the values of the interference terms aJ^T^^^ 

and crj/^'^^, the stronger the constraint on crj The fact that in the present data the 

measured asymmetries A lt and Att are consistent with zero implies that the interference 
terms crj^^^^ and a'^^^^^ are also consistent with zero. In this case, the inequalities (10) 
are trivially satisfied and give no information on o"2 



10.3 dependence of da^*P^^P/dMx 

Figure 7 shows the cross-section da"'*^^^^/ dMx as a function of for different Mx and 
W values. The data are also presented in Tables 3 and 4. The present measurements 
are shown together with the previous ZEUS results at low [7] and high [9]; the latter 
have been corrected for the residual double-dissociative background, taken to be 31%, as 
determined [9] by comparing those data with the LPS results [11] ^. The present results are 
consistent with the earlier ZEUS measurements [7, 9] and cover a wider kinematic region; 
notably, they reach higher values of Mx, lower values of Q^, as well as values of W close to 
the kinematic limit. The points at Mx = 5 GeV are all in the diffractive-peak region, since 
xjp < 0.01. The other bins have contributions from a;jp > 0.01. In all regions of X]p, the 
data exhibit a behaviour qualitatively similar to that of the total 7*p cross section, £7^/': 
da'^*P~'^P/dMx falls rapidly with at high Q^; as 0, the cross-section dependence 

on becomes weak, with da"'*'^^'^'^ / dMx approaching a constant, as expected from the 

^ As discussed in Section 10.5, a higher fraction was measured at lower [7]. 
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conservation of the electromagnetic current. The behaviour of da'^ ^^-^^ / dMx cannot be 
fitted with a simple, form-factor-like function of the type 1/ (Q^ + M].)"^, but is described 
by a pQCD-based model [48] at xjp < 0.01 and large (see the continuous curves in 
Fig. 7). This model and the comparison between its predictions and the data, as well as 
the curves on Fig. 7, are discussed in Section 10.7. 

A direct comparison between da'^*^~*^^/ dMx and the total photon-proton cross section is 
shown in Fig. 8, where the ratio {Mxda^*^""^^ / dM"^) /a^^^ is presented as a function of 
at different Mx and W values. The values of the 7*p total cross section were obtained 
from the ALLM97 parameterisation [49], which is consistent with the latest HI and ZEUS 
F2 data [50]. The plot shows that, in spite of their qualitative similarity, da'^ ^^^^/dMx 
and the total cross section exhibit some differences in their dependences. In the bin 
at Mx = 5 GeV, which has data from the diffractive peak region only, the ratio grows 
slowly with for < M"^ and then falls in the region dominated by qq fluctuations of 
the photon. At higher Mx values, which correspond to xjp larger than 0.01, is always 
smaller than Mj- and the ratio grows, indicating a softer dependence of da'^ v^^v j (IM\ 
than of o"7o/'- Here again, the low-xjp data in the high-Q^ region can be described by 
pQCD-based models of diffraction, as argued in Section 10.7, where the curves on Fig. 8 
are discussed. 



10.4 The structure function 

The data of Fig. 7 are presented in Fig. 9 (see also Tables 5-9) in terms of the structure 
function F^^^\ evaluated under the assumption that = 0. As discussed in Section 2, 
F-^^^^ is defined in this paper as the integral of F-^^^^ over the range < |t| < 1 GeV^. 
The figure shows X]pF2^^\(3, Q^, xjp) as a function of xjp for different values of (3 and Q"^. 
In order to maximise the kinematic overlap between the low-Q^ and the high-Q^ samples, 
only a subset of the low-Q^ sample is presented in Fig. 9 and Table 5. The values of 
X]pF^^^^ decrease with xjp a± small xjp, indicating that ^2^''^'' falls with xjp faster than 
1/xjp. At larger xjp, X]pF^^^^ flattens and, in some bins, increases with xjp, indicating 
a softer dependence of F^^^^ on xp. The xjp dependence can be parameterised in terms 
of the Pomeron intercept ajp(O). The extraction of «jp(0) is presented in Section 10.6, 
where the curves in Fig. 9 are discussed. 

The dependences of the structure function on and /? are presented in Figs. 10-12 
for different values of xjp. The structure function rises with in all of the explored 
kinematic region. Figure 10 compares the low-Q^ and high-Q^ results in the xp-jS region 
where they overlap. Between the low- and high-Q^ data, i.e. between (Q^) ~ 0.15 GeV^ 
and (Q^) ~ 10.5 GeV^, the increase is about a factor five; this steep rise, F^ oc Q^, reflects 
the flattening of da'^*P^^P/dMx for ^ (see Eqs. (1), (4)), and is a consequence of 
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the conservation of the electromagnetic current. The dependence becomes slower in 
the high-Q^ data (Fig. 11). In this region, the rise can be interpreted as a manifestation 
of QCD evolution; these positive scaling violations are due to the large gluon contribution 
to F2 (see Section 10.7.3). The behaviour is similar in all xjp bins, both in the diffractive- 
peak region and at larger xjp values, suggesting that the QCD evolution of the diffractive 
PDFs is largely independent of xjp. The (3 dependence of Fg^'"^'' (Fig. 12) has instead a 
different behaviour at different values of xp: at small xjp, ¥2^^^ has a weak (3 dependence 
with a tendency to rise at large values of /3; for xjp > 0.02, ¥2'^'^^ decreases as (3 approaches 
unity, as expected in a hadron (bearing in mind that (3 is the equivalent of Bjorken-x for 
the exchange). The solid curves shown on Figs. 11 and 12 are the predictions from a 
pQCD-based model [51] valid in the diffractive region and discussed in Section 10.7.2. 

The structure function Fg^*-^'* and the inclusive proton structure function F2 are compared 
in terms of the ratio X]pF^^^\x]p, x, Q'^)/F2{x, Q"^) calculated at fixed values of xjp. The 
values of F2(x, Q^) were obtained from the ALLM97 parameterisation. The ratio is pre- 
sented in Figs. 13 and 14 as a function of and x, respectively. The x range covered 
by the data is 5 x 10^^ < x < 5 x 10^^. The ratio is largely Q^-independent, with possi- 
bly some structure at high (3 and low xjp. A Q^-independent ratio would indicate equal 
scaling violations in the proton for the reaction 7*p — Xp and for inclusive DIS. The 
ratio grows with x in the diffractive-peak region, suggesting that the x dependence of the 
proton PDFs is different when the proton is probed in 7*p —>■ Xp and in inclusive DIS. 
At higher xjp, xp > 0.02, the ratio becomes flatter. 

10.5 Comparison with the results of the Mx method and of HI 

In this section, the present results are compared with those obtained with the Mx 
method [9] and with the HI results obtained with the rapidity-gap technique [14]. 

The Mx-method analysis [9] includes events in which the proton diffractively dissociates 
into a system N of mass M^v ^ 5.5 GeV. In order to facilitate the comparison, the present 
data were replotted using the binning of the M^-method publication [9]. 

The Mx-method points are higher than those obtained with the LPS method. The differ- 
ence was quantified by means of BEKW-type fits to the two data sets (see Section 10.7.1), 
which give a ratio of the M^-method to LPS points of Rmx = 1.55 ± O.OSlo;^^, where 
the first error includes the statistical and systematic uncertainties and the second is due 
to the normalisation uncertainty. Figure 15 shows xjpF-^^^^ for the two data sets after 
scaling down the Mx-method points by Rm^. The agreement between the LPS points 
and the renormalised Mx-method results is good, indicating that the difference is mainly 
in the normalisation. The normalisation difference can be attributed to the residual p- 
dissociative background in the Mx method. This background was estimated in previous 
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studies [9] by comparing the Mx results [9] and the earher LPS data [11]; there the ratio 
Rmx was found to be Rm^ = 1.451q;23, consistent with the present result. A similar study 
was performed [7] for the BPC region, resulting in Rmx = 1.85 ±0.38 (stat.), which is also 
consistent with the present measurement. The measured value of Rmx corresponds to a 
percentage of p-dissociative events in the sample of -Rdiss = {^ — ^/Rmx) = [35.5±3.3l7;^]%; 
the first error corresponds to the statistical and systematic uncertainties in quadrature, 
the second is due to the normalisation uncertainty. 

The agreement between the LPS and the Mx-method results, after taking the proton- 
dissociative background into account, lends support to the assumptions on which the Mx 
method is based. 

Figure 16 shows a comparison of the present F^^^^ results with those of the HI collabo- 
ration obtained with the rapidity-gap selection [14]. These data include a p-dissociative 
contribution with M^r < 1-6 GeV. The data are plotted in terms of xjpF^^^^ as a function 
oi xjp in different (3 and bins. The ZEUS points were extrapolated to the HI bin 
centres using the measured dependences. At small X]p, xjp < 0.01, the agreement is good, 
although with a tendency for the HI points to be higher than the present results at high 
Q^. While a normalisation difference is consistent with the presence of a p-dissociative 
contribution in the HI data, a dependence of this difference is not expected. The 
comparison indicates that the HI data have a stronger dependence than the present 
data. For xjp > 0.01, the HI data are also higher, but the shape is somewhat different, 
with a larger Reggeon-like contribution in the HI data. 

10.6 Extraction of the Pomeron intercept 

In the framework of Regge phenomenology, the xjp dependence of F^^^^ is related to the 
intercept of the Pomeron trajectory, the parameter that drives the energy dependence of 
the total hadron-hadron cross section at high energies [3]. The Pomeron intercept has 
been determined to be imGitomll 

[52] in soft hadronic interactions. However, the same 
parameter is significantly larger in the diffractive production of heavy vector mesons, 
notably in J/ip photoproduction (see e.g. [2]), reflecting the rapid rise of the cross section 
with W. This is a consequence of the increase of the parton densities in the proton 
at low X, which drives the rise of the cross section with decreasing x, and hence with 
decreasing xp (since xjp oc oc x). It is interesting to determine if such a deviation 

from the behaviour of the hadron-hadron data is also apparent in the inclusive diffractive 
dissociation of virtual photons. 

The high-Q^ data of Fig. 9 with xjp < 0.01 were fitted to the form: 

F,''^'\p,Q',xjp) = Mxjp) ■ Ff{(3,Q'), 
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i.e. assuming "Regge factorisation" (see Section 7). The Pomeron flux was parameterised 



The parameter a'jp was set to 0.25 GeV~^, consistent with the hadron-hadron data [53]. 
The parameter was taken to be 4.67 GeV~^, such that the relation h = h^—2a']p In xp [3] 
reproduces the results of Section 10.1. The values of F:f{P,Q^) in each /3 and bin 
and the Pomeron intercept (assumed to be (3 and independent) were treated as free 
parameters. The resulting Pomeron intercept is 



Varying a'jp in the range < a'jp < 0.4 GeV~^ causes ap^O) to change by laoso! chosen 
range includes the small values of a'jp measured in the diffractive photopro duct ion of J/i/j 
mesons [54]. A variation of between and 1 produces a change of aip(O) of +0.017. 
The quahty of the fit is good, with x^/ndf = 14.2/22 (considering statistical uncertainties 
only). The fact that the same value of ajp(O) fits the whole Q^, P region covered by the 
data indicates that, within the present accuracy, the hypothesis of Regge factorisation is a 
good approximation. The result does not change if the fit is extended to the low-Q^ region; 
in this case x^/nc?/ = 14.6/28 (statistical uncertainties only). The result of this latter fit 
is shown in Fig. 9. The extrapolation of the fit for xjp > 0.01, where the contribution 
from the exchange of the Reggeon trajectory becomes important, is also shown; indeed 
the fit does not describe this region satisfactorily. The present value of ajp(O) is consistent 
with that from HI, ajp(O) = 1.203 ± 0.020(stat.) ± 0.013(syst.)+|!;j35(model) (measured 
in the region 4.5 < < 75 GeV^) [14], and with the earlier ZEUS result, ttjp(O) = 
1.157 ± 0.009(stat.)+S;Si2(syst.) (measured in the region 7 < < 140 GeV^) [9]; it is 
higher than that of the soft Pomeron, suggesting that the parton densities probed in 
inclusive diffractive ep interactions also increase rapidly at small x and that a single 
Pomeron trajectory cannot simultaneously describe the high-Q^ diffractive data and the 
soft hadron-hadron data. 

10.7 Comparison with models 

As discussed earlier, the diffractive dissociation of virtual photons can be described in 
pQCD since the virtuality of the photon provides a hard scale. In the proton rest frame, 
the reaction can be viewed as the sequence of the photon fluctuating into a qq (or qqg) 
colour dipole, the dipole scattering off the proton and producing the final state X. At 
high centre-of-mass energies, these processes are widely separated in time. The qq, qqg 
fluctuations are described in terms of the photon wavefunction derived from QCD. The 



as [1] 




ajp(O) = 1.16 ± 0.02 (stat.) ± 0.02(syst.). 
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interaction of the dipole with the proton is mediated, in the lowest order, by the exchange 
of two gluons in a colour- singlet state. 

Several models of inclusive diffraction are available, which are discussed in review arti- 
cles [2,3]. The discussion of this section is restricted to two approaches based on the 
framework just outlined. The data were fitted (Section 10.7.1) with a parameterisation 
based on the model of Bartels et al. (BEKW) [48], which gives a satisfactory description 
of the earlier ZEUS results [7,9]. It is interesting to see if the same parameterisation is 
able to describe the present data which cover a wider kinematic region; the fit is also a 
useful tool to compare the present data and those based on the ZEUS analysis using the 
Mx method [9]. 

In Section 10.7.2, the results of this paper are also compared with the Golec-Biernat and 
Wiisthoff model based on the idea of the saturation of the dipole-proton cross section [55], 
which successfully describes both the inclusive ep scattering data and earlier diffractive 
data. 

Finally, Section 10.7.3 describes the results of an NLO QCD fit to the present high-Q^ 
data. In this approach, the dependence of the data is interpreted as due to the QCD 
evolution of the diffractive PDFs. A parameterisation of the diffractive PDFs at a starting 
scale is evolved according to the QCD evolution equations and fitted to the data. 

10.7.1 BEKW fit 

In the BEKW model [48], the dominant (leading-twist) contributions to the diffractive 
cross section in the kinematic domain of the present measurement come from fluctuations 
of transversely polarised virtual photons into either qq or qqg states. The (3 (and hence 
Mx) spectra of these two components are determined by general properties of the photon 
wave-function, with the qq contribution to the cross section proportional to (3{1 — (3) and 
the qqg contribution proportional to (1 — (3)^ , where 7 is a free parameter. For small 
values of Mx, the qq states dominate, while at large masses the qqg contribution becomes 
dominant. The model does not fix the xp dependence of the qq and qqg contributions, 
but assumes for both a power-like behaviour, x^^^ \ where the exponent n is determined 
from fits to the data. More explicitly, in the BEKW approach, the diffractive structure 
function can be parameterised as 

xjpF?'^^\(3, XTP, Q2) = ctF];^ + clF^^ + c,Fj^, 

where 

F,^-=(xo/x^)"-(«')-/5(l-/3), 
= (xo/x^)"«('3') ■ In (1 + Q'/Ql) ■ (1 - py. 
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The contribution of longitudinal photons, F^^, which is relevant only at high (3, was 
neglected in this analysis. A higher-twist term for qq states produced by transverse 
photons was also neglected. In the original BEKW model [48], the exponents fiT^giQ"^) 
are parameterised as riT^giQ"^) = n^'^ + n^'^ln (1 + InQ"^ /Ql)] here, this dependence was 
modified to riT^giQ'^) = riQ'^ + n^'^ In (1 + Q^/Ql), which is well defined also when 0. 

A fit was performed to the present high-Q^ data using the parameterisation described 
above. The fit was limited to the region xjp < 0.01, well within the diffractive peak. The 
parameters Ql, xq and tiq'^ were taken to be 0.4 GeV^, 0.01, and 0.13, respectively ^. 
The coefficients ct, Cg, ri['^ and 7 were determined in the fit, and have the following 
values: ct = 0.072 ± 0.006(stat.), Cg = 0.008 ± 0.001(stat.), n]^'^ = 0.053 ± 0.014(stat.), 
7 = 12.78 ± 2.08(stat.). The main features of the data are broadly reproduced by the 
fit, as shown in Fig. 7; the description of the dependence of the diffractive to the 
inclusive cross-section ratio is also reasonable, as seen in Fig. 8. This indicates that the 
framework in which the incoming virtual photon fluctuates into a colour dipole is, in 
general, adequate to describe diffractive processes in ep collisions. At the same time, the 
data suggest the increasing importance of the contribution from qqg states at low Q^, 
as indicated in Fig. 7. The fit gives only a qualitative description of the low-Q^ sample, 
which is outside the region of applicability of pQCD; these points were not included in the 
fit. The fit is also lower than the high-Q^ data in the high-Mx bins that have xp > 0.01- 
0.02, suggesting that different mechanisms, such as Reggeon exchange, are at work in the 
diffractive-peak region and at high xjp. In this region the discrepancy between the data 
and the fit can be taken as an estimate of the contribution to the cross section due to 
exchanges other than the Pomeron. 

10.7.2 Saturation model 

In the saturation model by Golec-Biernat and Wiisthoff [55], diffractive DIS is also de- 
scribed in terms of the interaction of the qq (qqg) fluctuation of the virtual photon with 
the proton. At high Q^, the dipole-proton cross section is obtained from pQCD and is 
proportional to the square of the transverse size of the dipole, which is in turn propor- 
tional to As decreases, the rise of the cross section with the dipole size would 
violate unitarity and is tamed by requiring that it saturates at a typical value of the 
hadron-hadron cross section. The value of at which saturation occurs is x-dependent. 
The parameters of the model were obtained from a fit to F2 data. The latest modification 

5 Since titJQ'^) = r^^^'^ + nf'''ln(l + QVQg), uq gives the xp dependence of F^'^^^ for Q'^ 0. If 
F2^'^^ is assumed and aip{t) is taken to be the soft Pomeron trajectory (aip(O) = 1.096, 

a'p = 0.25 GcY-'^), then no = 0.13. 
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of the model [51], denoted by BGK in the following, includes the QCD evolution of the 
gluon distribution. 

Figures 11, 12 and 17 show the comparison of the measured structure function F^^^'' with 
the BGK prediction [51]. In the region of apphcability of the model, xjp < 0.01 and 
larger than a few GeV^, the xjp, and (3 dependences of F^^^^ are adequately described, 
although the data are slightly higher than the model prediction. The extrapolation of the 
model to large Xjp values, beyond the Pomeron-dominated region, is significantly lower 
than the data, and the discrepancy increases with xjp; in this region, the (3 dependences 
of the data and the model are also markedly different. 

Both the BEKW and the saturation model imply Regge factorisation breaking. The fact 
that these models describe the data is not in contradiction with the possibility to fit the 
same data assuming Regge factorisation, as was done in Section 10.6, since the magnitude 
of the predicted violation is smaller than the precision of the present data. 



10.7.3 QCD fit 

An NLO QCD fit was performed to the present high-Q^ data together with the recent 
ZEUS results on diffractive charm production in DIS [56]. The latter are important to 
constrain the gluon contribution to the diffractive PDFs. The fit was limited to the 
data in the region > 2 GeV^ and xjp < 0.01. Regge factorisation was assumed, 
F^^^\(3,Q'^,X]p) = fip{xip) ■ F^{(3,Q'^) (see Section 7), and the Pomeron flux was taken 
to be of the Donnachie-Landshoff form [53] 

where /5o = 1.8 GeV"\ Fi{t) is the elastic form factor of the proton, ajp(O) was fixed to 
the result given in Section 10.6 and a'p was set to 0.25 GeV~^. The results do not change 
if Eq. (9) is used. 

The diffractive parton distributions (quark flavour singlet and gluon) were parameterised 
at the starting scale, Ql = 2 GeV^, using the general polynomial form zf\z) = (ai + 
a2Z + asz^) ■ (1 — z)""^, where z is the parton fractional momentum. For the light quark 
distribution, it was assumed that u = d = s = u = d = s; it was verified that setting 
the strange quark density to zero at the starting scale produces no appreciable change in 
the results. Charm quarks were treated in the Thorne-Roberts variable flavour number 
(TRVFN) scheme [57], with the charm-quark mass, m^, set to 1.45 GeV. The NLO evolu- 
tion package QCDNUM [58] was used to evolve the PDFs from the starting scale to the 
values of each data point. The evolved PDFs were then fitted to the data. 
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The result of the fit is shown by the lines in Fig. 18. They satisfactorily reproduce the 
measurements, with /ndf = 37.8/36 (statistical errors only). The resulting fraction of 
the t-channel momentum carried by gluons is (82 ± 8(stat.)]'i^g(syst.))% at = 2 GeV^, 
consistent with earlier ZEUS [59] and HI [14] results, but higher than that found in a 
recent QCD analysis of the same data by Martin, Ryskin and Watt [60]. The system- 
atic uncertainty includes the contributions listed in Section 9; in addition, the charm 
quark mass was varied between 1.3 and 1.6 GeV, and the relative normalisation between 
the charm and the F^^^^ data was changed by li3%, reflecting the uncertainty on the 
proton-dissociative background and on the luminosity in the charm data, as well as the 
normalisation uncertainty of the LPS data. The fixed-fiavour-number scheme (FFNS) 
was used instead of TRVFN, without any significant change of the results. Various PDF 
parameterisations at the starting scale were tried, including the function used by HI [14]. 
The shape of the fitted PDFs changes significantly depending on the functional form of 
the initial parameterisation, a consequence of the relatively large statistical uncertainties 
of the present sample. Therefore, these data cannot constrain the shapes of the PDFs. 
However, the integrals over z of the fitted PDFs and notably the fraction of the t-channel 
momentum carried by gluons are robust and change only slightly with the parameterisa- 
tion chosen. This contribution was included in the systematic uncertainties quoted. 

Also shown in Fig. 18 is the ratio of the charm contribution to the diffractive structure 
function, p^^^^'^^ [56]^ and the present F^^^^ results as a function of /? (see Table 10). 
The ratio increases with increasing and decreasing (3, up to values of 30%. The ratio 
is well described by the fit. 



11 Summary 

New measurements have been presented of the reaction e~^p — >■ e~^Xp in the regions 
0.03 < < 0.60 GeV^ and 2 < < 100 GeV^ The scattered proton was measured 
in the ZEUS leading proton spectrometer, and was required to carry a fraction xl of the 
incoming proton momentum of at least 90%. The data cover the region 0.075 < \t\ < 
0.35 GeV2. 

The results can be summarised as follows: 

• the t dependence of the cross section is exponential, with a t-slope b = 7.9±0.5(stat.)]|iQ5(syst.) 
GeV~^ for xp < 0.01. The slope is independent of but decreases with xjp; 

• there is no observed $ dependence of the cross section, indicating that the interfer- 
ence terms between the longitudinal and transverse amplitudes and between the two 
transverse amplitudes are consistent with zero in the measured kinematic region; 
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• the cross-section da^'P^^^/dMx falls rapidly with at high but approaches a 
constant as ^0. This behaviour is similar to that of the total photon-proton cross 
section, and is a consequence of the conservation of the electromagnetic current. In 
detail, it was found that the cross-section da"'*^^'^^ / dMx falls with more slowly 
than a2ot when < M]^ and faster than cr^t^ when Q"^ > M|; 

• the data were also analysed in terms of the structure function F^^^\ The xjp, and 
f3 dependences of F^^^^ were studied; 

— F^^^^ falls with xjp faster than 1/xjp for Xjp < 0.01 and more slowly at larger 
values of xjp. From the xjp dependence of F^^^^ at low xjp, the Pomeron intercept 
q;jp(0) was measured to be ajp(O) = 1.16 ± 0.02 (stat.) ± 0.02(syst.), higher than 
that of the soft Pomeron, and similar to that measured in the photoproduction of 
heavy vector mesons. This suggests that, also in the present reaction, the virtual 
photon probes the proton in a region where the parton density increases quickly 
with decreasing x. 

— F^^^^ rises with over the whole measured region. The increase is very signif- 
icant, about a factor five, between the low-Q^ and the high-Q^ region. In the 
high-Q^ region, the rise becomes softer, and is reminiscent of the logarithmic scal- 
ing violations of the proton structure function. Positive scaling violations reflect a 
large gluon density. This is confirmed by an NLO QCD analysis of the present data 
for xjp < 0.01 in conjunction with the earlier ZEUS results on diffractive charm 
production [56]. The analysis indicates that the fraction of the t-channel momen- 
tum carried by gluons is (82 ± 8(stat.)]'^^g(syst.))% at = 2 GeV^. However, the 
present data are not precise enough to constrain the shapes of the PDFs. 

— The (3 dependence of F^^^^ changes with xp. For xp < 0.01, F2^^^ grows with 
(3. For values of xjp > 0.01-0.02, Fj^''^'* decreases with (3. The latter behaviour 
is similar to that of the structure functions of hadrons as a function of x, and is 
consistent with the hypothesis that, at large xjp, the t-channel exchange mediating 
the photon-proton interaction is a meson-like object; 

• the results presented are consistent, in the small xp region, with the predictions of 
pQCD-based models of diffraction. In particular, the data were compared with models 
in which the virtual photon fluctuates into qq or qqg colour dipoles which then interact 
with the proton via the exchange of a gluon pair. 
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0.06 


4.19 ± 0.40+°;^^ 



Table 1: Fitted values of the t-slopes. The first uncertainty is statistical, the 
second systematic. 
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Table 2: Fitted values of the asymmetry parameter Alt- The first uncertainty is 
statistical, the second systematic. 
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Table 3: The diffractive cross section, da'^ ^^^^/dMx, for the low-Q"^ sample. 
The first uncertainty given is statistical, the second systematic. 
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Table 4: The diffractive cross section, da'^'^^^^/dAdx, for the high-Q^ sample. 
The first uncertainty given is statistical, the second systematic. 
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Table 5: The dijfractive structure function multiplied byxjp, xjpFg Q^, xjp), 
for the low-Q'^ sample. The first uncertainty given is statistical, the second system- 
atic. 
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Table 6: The diffractive structure function multiplied by xjp, X]pF2 {(3, Q^, xjp), 
for the low-Q^ sample. The first uncertainty given is statistical, the second sys- 
tematic. These are the points plotted in Fig. 12. The data in this table are not 
independent of those in Table 5. 
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Table 7: The dijfractive structure function multiplied hyxjp, xjpFg {(3^Q^,X]p), 
for the high-Q^ sample, part I. The first uncertainty given is statistical, the second 
systematic. 
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Table 8: The diffractive structure function multiplied by xjp, xjpF^ {(3, Q^, xjp), 
for the high-Q"^ sample, part II. The first uncertainty given is statistical, the second 
systematic. „„ 
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Table 9: The diffractive structure function multiplied by xjp, xjpFg (/?, Q^, xp), 
for the high-Q"^ sample, part III. The first uncertainty given is statistical, the second 
systematic. 
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Table 10: The ratio of F2 [56] and the present F2 measurement. The 
first uncertainty given is statistical, the second systematic. 
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Figure 2: Comparison of the measured (points) and Monte-Carlo simulated (his- 
tograms) distributions for Xl, \t\, , W, Mx and xp in the low-Q"^ analysis. The 
and W distributions were obtained without the LPS requirement (see text). 
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Figure 3: Comparison of the measured (points) and Monte-Carlo simulated (his- 
tograms) distributions for xl, \t\, Q^, W, Mx and xjp in the high-Q"^ analysis. 
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Figure 4: Distribution of E + pz for the high-Q'^ events. The hatched histogram 
represents the beam-halo sample obtained as discussed in the text. The empty his- 
togram is the sum of the RAPGAP Monte Carlo and the beam-halo contribution. 
The vertical dashed line is at E + pz = 1655 GeV, the value of the cut used to 
suppress beam-halo events. 
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Figure 5: (a) The differential cross-section da'^^"*'^^^ /dt in the region xp < 
0.01, 2 < < 100 GeV^ and Mx > 1.5 GeV. The inner error bars show the 
statistical uncertainties and the full bars indicate the statistical and the systematic 
uncertainties added in quadrature. The overall normalisation uncertainty of ±10% 
is not shown. The line shows the result of the fit described in the text, (b) The 
value of the slope parameter h of the differential cross-section da'^^^'^^^ / dt as a 
function of Q^. (c) The value of the slope parameter b of the differential cross- 
section da'^'^~"^^^ / dt as a function of xjp. The mean value of (3 in each bin is also 
given, (d) The value of the slope parameter 6p2 of the differential cross-section 
da'^P~*'^^P / dp"^ as a function of xp. The symbols labelled ZEUS 97 indicate the 
present results. Earlier ZEUS results are also shown: ZEUS 95 [8], ZEUS 94 [12] 
(Q^ = 0) and ZEUS 94 [11] (Q'^ >0). 
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Figure 6: The differential cross-section da'^^^'^^^ / in the kinematic ranges (a) 
0.00025 < xjp < 0.01 and (h) 0.01 < xp < 0.07. The error bars show the statistical 
uncertainty. The line shows the result of the fit described in the text. The azimuthal 
asymmetry Alt cls a function of (c) xjp, (d) (3, (e) \t\ and (f) Q"^ for xp < 0.01. 
The inner error bars show the statistical uncertainties and the full bars indicate the 
statistical and the systematic uncertainties added in quadrature. 
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Figure 7: The cross-section da"**^^-^^ / dMx as a function of Q'^ at different Mx 
and W values. The inner error bars show the statistical uncertainties and the full 
bars indicate the statistical and the systematic uncertainties added in quadrature; in 
several bins the size of the bars is smaller than that of the symbol used for the points. 
The overall normalisation uncertainty of tiQ% is not shown. The symbols labelled 
ZEUS 97 indicate the present results. Earlier ZEUS results are also shown: ZEUS 
96-97 Mx [7j, ZEUS 94 Mx [9]. The solid lines are the result of the BEKW fit to 
the present high-Q^ data, described in Section 10.7.1; the dashed lines indicate the 
extrapolation outside the fit region. The dotted lines indicate the qqg contribution. 
The xp ranges given refer to the coverage of the data. 
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Figure 8: The ratio of Mjrda'^ ^""^^/dMjr to the total virtual photon proton 
cross section as a function of at different Mx and W values. The inner error 
bars show the statistical uncertainties and the full bars indicate the statistical and 
the systematic uncertainties added in quadrature. The overall normalisation un- 
certainty of~^\l% is not shown. Earlier ZEUS results are also shown: ZEUS 96-97 
Mx [7], ZEUS 94 Mx [9]. The solid lines are the result of the BEKW fit divided 
by the ALLM97 parameterisation, as described in Section 10.7.1. The dashed lines 
indicate the extrapolation outside the fit region. 
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Figure 9: The diffractive structure function multiplied by xjp, X]pF^^^^ Q^, xjp), 
as a function ofxp, for different values of (3 and . The inner error bars show the 
statistical uncertainties and the full bars indicate the statistical and the systematic 
uncertainties added in quadrature. In some bins the size of the bars is smaller 
than that of the symbol used for the points. The overall normalisation uncertainty 
of ^10% is not shown. The vertical dotted lines indicate xp = 0.01. For the 
low-Q"^ points, the average value of Q"^ in each fi-xp bin varies between 0.10 and 
0.22 GeV^ (see Table 5). The solid lines show the result of the Regge fit described 
in Section 10.6. The dashed curves are the extension of the fit for xjp > 0.01. 
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Figure 10: The dijfractive structure function multiplied by X]p, 

xipF^^^\l3,Q'^,xip), as a function of , for the low- and high-Q"^ data, in 
the (3-xp region in which the two data-sets overlap. The inner error bars show the 
statistical uncertainties and the full bars indicate the statistical and the systematic 
uncertainties added in quadrature. The overall normalisation uncertainty of 
is not shown. The data shown here are those presented in Table 5 rebinned in xl 
slightly for plotting purposes. 
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Figure 11: 

,D(3) 



The diffractive structure function multiplied by Xjp, 
xipF^^'^'{/3,Q'^,X]p), as a function of Q^, for different values of xjp and (3. 
The inner error bars show the statistical uncertainties and the full bars indicate 
the statistical and the systematic uncertainties added in quadrature. The overall 
normalisation uncertainty o/1Jq% is not shown. The solid lines are the prediction 
of the saturation model of Bartels et al. [51] (BGK) discussed in Section 10.7.2; 
the dashed lines indicate the extrapolation of the model beyond Xjp = 0.01. 
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Figure 12: The diffractive structure function multiplied by Xp, 

xipF^^^\(3,Q'^,X]p), as a function of (3, for different values of xjp and Q^. 
The inner error bars show the statistical uncertainties and the full bars indicate 
the statistical and the systematic uncertainties added in quadrature. The overall 
normalisation uncertainty of tiQ% is not shown. The low-Q^ data in the bin 
labelled xjp = 0.06 have xjp = 0.067 and xjp = 0.066, respectively (see Table 6). 
The solid lines are the prediction of the saturation model of Bartels et al. [51] 
(BGK) discussed in Section 10.7.2; the dashed lines indicate the extrapolation of 
the model beyond xp = 0.01. 
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Figure 13: The ratio of the diffractive to the inclusive structure functions, 
xipF^^^\x]p, x,Q'^)/ F2{x,Q^) , as a function of at different values of xp and 
X. The values of F2{x,Q^) were obtained from the ALLM97 parameterisation. The 
inner error bars show the statistical uncertainties and the full bars indicate the 
statistical and the systematic uncertainties added in quadrature. The overall nor- 
malisation uncertainty of ^j;o% is not shown. The horizontal lines indicate the 
average value of the ratio and are only meant to guide the eye. 
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Figure 14: The ratio of the diffractive to the inclusive structure functions, 
X]pF^^^\xip,x,Q'^)/F2{x,Q'^), as a function of x at different values of xp and 
. The values of F2{x,Q^) were obtained from the ALLM97 parameterisation. 
The inner error bars show the statistical uncertainties and the full bars indicate 
the statistical and the systematic uncertainties added in quadrature. The overall 
normalisation uncertainty c»/lio% ^■^ shown. The horizontal lines indicate the 
average value of the ratio and are only meant to guide the eye. 
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Figure 15: The diffractive structure function multiplied by xjp, 

X]pF^^^\l3,Q'^,X]p) as a function of xjp, for different values of (3 and 
for the LPS and Mx -method [9] analyses; the latter points are rescaled by 
1/Rmx = 0.645, as discussed in the text. The inner error bars show the statistical 
uncertainties and the full bars indicate the statistical and the systematic uncer- 
tainties added in quadrature. The overall normalisation uncertainties of 
(LPS data) and 2% (Mx-method data) are not shown. 
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Figure 16: The diffractive structure function multiplied by Xjp, 

X]pF^^^\(3,Q'^,xip) as a function of xjp, for different values of (3 and 
for the LPS and the HI data [14]- The inner error bars show the statistical 
uncertainties and the full bars are the statistical and the systematic uncertainties 
added in quadrature. The overall normalisation uncertainties oftio% (LPS data), 
±6% (HI data, < 9 GeV^) and ±4.8% (HI data, > 9 GeV^) are not shown. 
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Figure 17: 

,D(3) 



The diffractive structure function multiplied by Xjp, 
X]pF^^'^'{j3,Q'^,xip), as a function of xjp, for different values of (3 and Q'^ . 
The inner error bars show the statistical uncertainties and the full bars indicate 
the statistical and the systematic uncertainties added in quadrature. The overall 
normalisation uncertainty o/1Jq% is not shown. The solid lines are the prediction 
of the saturation model of Bartels et al. [51] (BGK) and are discussed in Sec- 
tion 10.7.2. The dashed curves are the extrapolation of the model prediction for 
Xjp > 0.01. 
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Figure 18: (a) The diff'ractive structure function multiplied by Xjp, 
xipF^^^\(3,Q^,xip), as a function of , for different values of xjp and j3. 
The inner error bars show the statistical uncertainties and the full bars are the 
statistical and the systematic uncertainties added in quadrature. The overall 
normalisation uncertainty of ^Jq^ is not shown, (b) The measured charm contri- 
bution to the diffractive structure function multiplied by xjp, X]pF^^^^'^^{j3,Q'^,X]p), 
as a function of (3, for different values of and xp = 0.004 [56]. (c) The ratio 
of F^^^^'^'^ [56] and the present F^^^^ measurement as a function of (3. The solid 
lines are the result of QCD NLO fit described in the text. 
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